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Abstract 

We use stellar structure dynamics arguments to extract bounds on the relevant pa- 
rameters of scalar field models: the putative scalar field mediator of a fifth force with 
a Yukawa potential, the new variable mass particle (VAMP) models, and a phenomeno- 
logically viable model to explain the Pioneer anomaly. We consider the polytropic gas 
model to estimate the effect of these models on the hydrostatic equilibrium equation and 
fundamental quantities such as the central temperature. The current bound on the solar 
luminosity is used to constrain the relevant parameters of each model. 

1 Introduction 

Scalar fields play a crucial role in particle physics and cosmology. Indeed, in inflation, the 
potential of a scalar field, the infiaton, acts as a dynamical vacuum energy that allows for an 
elegant solution of the initial conditions problems [1]. This prominent role of scalar fields is 
also evident in models to explain the late time accelerated expansion of the Universe in vacuum 
energy evolving and quintessence models [2], as well as in the Chaplygin gas dark energy-dark 
matter unification model [3]. Scalar fields have also been proposed as dark matter candidates 
[4]. Furthermore, it has been recently proposed that a scalar field can be also at the source of 
the anomalous acceleration detected by the Pioneer spacecraft [5]. 

Scalar fields may also have astrophysical implications as, for instance, the mediating boson of 
an hypothetical fifth force, which should yield the measurable effects on celestial bodies, besides 
the other known forces of nature. Although the origin of these fields is rather speculative, they 
can all be described by the Yukawa potential, written here as Vy(r) = Ae~"^^/r, where A is 
the coupling strength and m is the mass of the field, which sets the range of the interaction. 
Ay = m~^. 
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Models leading to a Yukawa- type potential can be found in widely distinct areas such as 
braneworld models, scalar-tensor theories of gravity and in the study of topological defects. In 
braneworld models, one considers our Universe as a 3-dimensional world-sheet embedded in a 
higher dimensional bulk space [6] . Symmetry considerations about the brane and its topological 
properties can be implemented to constrain the evolution of matter on the brane and gravity 
on the brane and in the bulk. 

Braneworld models are rather trendy in cosmology and allow, for instance, for a solution 
for the hierarchy problem, whether the typical mass scale of the bulk is comparable with the 
electroweak breaking scale. Mew ~ TeV. As a result, a tower of Kaluza-Klein (KK) massive 
tensorial perturbations to the metric appears. Following the KK dimensional reduction scheme, 
the masses (eigenvalues) of these gravitons (eigenfunctions) are ordered. Most braneworld 
models consider one first light mode with cosmological range, and hence all ensuing modes 
have sub-millimeter range. 

Given the relevance of scalar fields, the search for bounds on the Yukawa parameters is 
crucial, so to exclude unviable models and achieve some progress in the study on those that 
appear feasible. Most experimental tests of a "fifth" force have been conducted in the vacuum; 
in the authors opinion, a study on the way this force should affect stellar equilibrium is lacking, 
and constitutes one of the motivations for this work. 

The bounds on parameters A and Xy = include the following (see [7, 8] and references 
therein): laboratory experiments devised to measure deviations from the inverse-square law, 
sensitive to the range 10~^ m < Ay < 1 m, and constraining A to be smaller than 10~^; 
nucleosynthesis bounds which imply that A < Ax 10~^ for Ay < 1 m; gravimetric experiments, 
sensitive in the range of 10 m < Ay < 10^ m, suggesting A < 10~^; sateUite tests probing 
ranges of about 10^ m < Ay < 10^ m, showing that A < 10~^; and radiometric data of the 
Pioneer 10/11, Galileo and Ulysses spacecrafts suggesting the existence of a new force with 
parameters A = —10^'^ and Ay = 4 x 10^^ m [9]. It is striking that, for Ay < 10^'^ m and 
Ay > 10^'^ m, A is essentially unconstrained. Considerations on higher dimensional superstring 
motivated cosmological solutions hint that modifications to Newtonian gravity will occur in the 
short range region. Ay < 10~^ m (cf. Figure 2 below). This range also emerges if one assumes 
the observed vacuum energy density to be related with scalar or vector/tensor excitations [10]. 

In high energy physics, it is widely accepted that the mass of fermions results from the Higgs 
mechanism, in which a scalar field coupled to the right and left component of a particle acquires 
a vacuum expectation value that acts as a mass term in the Lagrangian density. This behaviour 
of the Higgs-scalar field depends on the presence of a potential which acquires non- vanishing 
minima and, therefore, cannot evolve monotonically. 
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Can one relax this last feature of the Higgs mechanism? This has been the main motivation 
behind the variable mass particle proposal [11]. In these models it is assumed that there are 
some yet unknown fermions whose mass results not from a coupling to the Higgs boson, but to a 
quintessence-type scalar field with a monotonically decreasing potential. This potential has no 
minima, yet the coupling of the scalar field to matter can be included in an effective potential 
of the following form V"e//(0) = V{(f)) + Xn^(f), where is the number density of fermionic 
VAMPs and A is their Yukawa coupling. In this way, a minimum is developed and the ensuing 
vacuum expectation value (vev) is responsible for the particle's mass. Since in a cosmological 
setting the density depends on the scalar factor a{t), this mass will vary on a cosmological time 
scale. 

Before proceeding, note that the present analogy is not perfect. Indeed, while the Higgs 
mechanism relies on a spontaneous symmetry breaking, where the vev experiences a transition 
from a vanishing to a finite value, the VAMP idea assumes that no vev exists if the matter 
term of the effective potential is "switched off", whereas it is always non- vanishing when the 
latter is considered. 

For definitiveness, we choose a potential of the quintessence-type form, V{(f)) = Uo4>~^, where 
Uq has dimensionality Mp~^^. The effective potential Ve//(0), acquires a vev given by 

/ ^ l/(l+p) 

Since the number density evolves as = n^o'2(^)~^i while (p evolves as = (l)oa^^^^^^\ 

where 4>q is the present value of the scalar field, then its mass is given by 



ml 



= [p{p + l)uocPo^^'^] a-3(2+^')/(^+^') , (2) 



and the mass of the VAMP fermions is [11] 

mv, = A0oa'/(i+^) . (3) 

Considering the evolution of the energy density contributions as a function of the redshift, 
z, one can compute the age of the Universe 

t= /^'^'"^ [l - fio + nuox-' + dx , (4) 

Jo a' Jo ^ J 

where Hq — 100 h Km Mpc"^, Qmo is the energy density of normal baryonic plus dark 
matter, Jly is the energy density due to the potential driving and is the total energy 
density of the Universe. The hmiting case where = — 1, ^mo — yields 
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to^lH^'{l+p-') . (5) 

Assuming that the VAMP particles, ijj, were relativistic when they decoupled from thermal 
equilibrium, then it follows that [11] 

my, = 12.7 nv,o/iV;V/2 el^ , (6) 

where is the ratio of Qeff, the effective number of degrees of freedom of ^z^, to g^f, the total 
effective number of relativistic degrees of freedom at freeze-out. Furthermore, in terms of the 
Yukawa coupling, 

1.0 = 1.02x10-9^1^61-5 , (7) 

and thus 

= 1.00 X lO-^^S-a^'^^ eV . (8) 

However interesting, VAMP models have not been subjected to a more concrete analysis 
mostly due to a significant caveat, namely, the introduction of exotic ip fermions, the VAMP 
particles, and the consequent derivation of the cosmologically relevant quantities in terms of 
their unknown relative density, Q^o- As a result of this somewhat arbitrary parameter, plus 
the unknown coupling constant with the scalar field and the potential strength, there is little 
one can do in order to draw definitive conclusions from VAMP models. 

The present study attempts to overcome this drawback, by asserting that, asides from the 
hypothetical existence of the assumed exotic particles, all fermions couple to the quintessence 
scalar field. In this assumption one considers that fermionic matter couples mainly to the Higgs 
boson, so that the VAMP mass term is a small correction to the mass acquired by the Higgs 
mechanism. Thus, exotic VAMP particles are defined by their lack of the Higgs couphng. 

The extension of the VAMP proposal to usual fermionic matter implies in a correction to 
the cosmologically relevant results obtained in Refs. [11, 12]. However, these corrections should 
be negligible, as one assumes that the "Higgs to quintessence" coupling ratio is small, that is, 
at cosmological scales the exotic VAMP particles dominate the VAMP sector of usual fermions. 

Notice that the vev resulting from the effective potential depends crucially on the particle 
number density n^. Hence, it is logical to expect that the effect of this variable mass term in 
a stellar environment should be more significant than in the vacuum. This sidesteps the model 
from the usual cosmological scenario with a temporal variation, to the astrophysical case with 
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an isotropic spatial dependence, thus allowing one to attain bounds on model parameters from 
known stellar physics observables. 

Another object of this study concerns the Pioneer anomaly. This consists in an anomalous 
acceleration inbound to the Sun and with a constant magnitude of qa — (8.5±1.3) x 10"^'^ ms~'^, 
revealed by the analysis of radiometric data from the Pioneer 10/11, Galileo and Ulysses space- 
craft. Extensive attempts to explain this phenomena as a result of poor accounting of thermal 
and mechanical effects and/or errors in the tracking algorithms were presented, but are now 
commonly accepted as unsuccessful [9]. 

The two Pioneer spacecraft follow approximate opposite hyperbolic trajectories away from 
the Solar System, while Galileo and Ulysses describe closed orbits. Given this and recalling 
that one has three geometrically distinct designs, an "engineering" solution for the anomaly 
seems not very plausible. 

Hence, although not entirely proven and even poorly understood (see e.g. Ref. [5] and 
references within), the Pioneer anomaly, if a real physical phenomenon, should be the manifes- 
tation of a new force. This force, in principle, acts upon the Sun itself, and thus lends itself to 
scrutiny under the scope of this study. It turns out that one can model its effect and constrain, 
even though poorly, the only parameter involved, the anomalous acceleration a^. 

2 The polytropic gas stellar model 

Realistic stellar models, arising from the assumptions of hydrostatic equilibrium and Newtonian 
gravity, rely on four differential equations, together with appropriate definitions [13, 14]. This 
intricate system requires heavy-duty numerical integration with complex code designs, and an 
analysis of the perturbations induced by a variable mass is beyond the scope of the present 
study. Instead, we focus the polytropic gas model for stellar structure: this assumes an equation 
of state of the form P — Xp^+^Z", where n is the so-called polytropic index, that defines 
intermediate cases between isothermic and adiabatic thermo dynamical processes, and K is 
the polytropic constant, defined below. This assumption leads to several scahng laws for the 
relevant thermodynamical quantities. 



p = , (a) (9) 

T^TM) , (6) 

P = PMT+' , (c) 

where pc, Tc and Pc are the values of the density, temperature and pressure at the center of the 
star. 
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The function 0, responsible for the scahng of P, p and T, is a dimensionless function of the 
dimensionless variable ^, related to the physical distance to the star's center hj r — a^, where 



a 



{n + l)K 
4nG 



(l-n)/n 
rc 



1/2 



(10) 



(11) 



and 



n + 1 



:{3-n)/n 



(n-l)/n)- 



(12) 



so that R is the star's radius, M its mass and defined by 9{^i) = 0, corresponds to the surface 
of the star (actually, this definition states that all quantities tend to zero as one approaches the 
surface. Its unperturbed value is ^f'^ — 6.89685, as given in Ref. [13]. 

The function 9{^) obeys a differential equation arising from the hydrostatic equilibrium 
condition 



the Lane-Emden equation: 



dr \dr p 



-a 



dM{r) 
dr 



(13) 



1^ ^2^\ ^ _nr, 



(14) 



We point out that the physical radius and mass of a star appears only in Eq. (11) so that the 
behaviour of the scaling function 9{^) is unaffected by their values. Hence, the stability of a star 
is independent of its size or mass, and different types of stars correspond to different polytropic 
indices n. This kind of scale-independent behavior is related to the homology symmetry of the 
Lane-Emden equation. 

The first solar model ever considered corresponds to a polytropic star with n — 2> and was 
studied by Eddington in 1926. Although somewhat incomplete, this simplified model gives rise 
to relevant constraints on the physical quantities. 
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3 Results 



3.1 Yukawa potential induced perturbation 

In this subsection we look at the hydrostatic equilibrium equation with a Yukawa potential: 



GM(r) [1 + Ae-"n , , , 
dP = -p{r)dr 



(15) 



which, after a small algebraic manipulation, implies that 



dr \dr p J 



l + Ae 



+ 



GM{r)Ame 



(16) 



The last term is a perturbation to the usual Lane-Emden equation, obtained by substituting 
r — and p — PcO"; one gets 



edA^ d^ - ^ 



+ 



M(^)ylme-"'"^ 



;i7) 



Since a = (^^^^pi'"")/'*)'^' and M{i) = -47r 
tion term can be written as 



{n+l)K 
4nG 



3/2 



p{3 n)/2n^2<m^ ^j^g secoud perturba- 



Anpca'^^^ 



Furthermore, one has 



{n + l)K 
AttG 



1/2 



(18) 



K = 



Pr = 



WnGM^/R'^ , 

i^p(n+l)/n ^ 



(19) 



where i?, M is the radius and mass of the star and iV„, Wn are numbers which depend on n, 
and for which one takes the tabulated values, valid for the unperturbed equation. Hence, 



Pc 



(20) 



Substituting into Eq.(18), one obtains 



{n-^l)K 
AttG 



1/2 



^(l-n)/2n^^^g-ame ^ _ ^!!±l7V;^/(n+l)|y a-n)/2(n+l)^^^^g-ame . (21) 



If one now defines the dimensionless quantities 
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C„= (^!^y^'Ar;^/("+i)iy^i-)/2("+i) , ^ = mR , (22) 
the perturbed Lane-Emden equation acquires the form 



1 d de 



-9' 



de 



ed(^ - ,^ + Ae-^\-rAC^-e-'^ . (23) 

One can ehminate a on the exponential term by writing —am^ — —amR^/R — —^C/Cs, 
where corresponds to the value at the surface of the star. Since one must specify it prior to 
integration of the differential equation, one assumes that ~ the latter being the tabulated 
value for ^(.^i) = 0. This is in good approximation, since the solar temperature is very small 
when compared at the surface to its central value, = 5.778 x 10^ K = 3.7 x 10~'^Tc. Hence, 
one gets 



1 d de 



de 



edidi - ,l + ^e-^«']-7^C,.-.-«. . (24) 

This notation makes clear that the perturbation vanishes for A ^ 0. Notice that, if m ~ 
R~^, then 7 ~ 1. However, since m should arise from a fundamental theory, this would be the 
case only for stars of a particular size = Ay. Since 7 = mR = R/Xy, it is clear that large 
stars (i? ^ Ay — >■ 7 ^ 1) are perturbed only within a small central region, where ^/^i -C 1. 

It is also apparent that the Yukawa-type perturbation breaks the invariance under homol- 
ogous transformations, since 7 = mR explicitly depends on R, the radius of the star. Hence, a 
bound on 7 obtained from the luminosity or other observables is, for a given m, equivalent to 
a bound on the maximum size of a polytropic star of index n. 

The boundary conditions are unaffected by the perturbation: from the definition p — Pce{^), 
one gets ^(0) = 1; the hydrostatic equation Eq. (16), in the hmit ^ — > 0, still imposes that 

M^Kl5=o = 0. 

By specifying A and 7, we may solve Eq. (16) numerically and get relative deviations from 
the unperturbed solution. Since a star's central temperature is given as a function of M and 
R and in terms of the mean molecular weight /i, the Hydrogen mass H and the Boltzmann 
constant k hy Tc — YnjiGM/ R, with 



■3|, 



l/n 

(25) 

one can compute the relative changes on Tc for different values of A and m. The results are 
presented in Figure 1. The parameters were chosen so that the Yukawa interaction Ay ranges 
from O.IR to lOR: in the first case, the interaction is mainly located in the interior of the star. 
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0.1 



Figure 1: Relative deviation from unperturbed central temperature T^/Tco — 1, for A ranging 
from 10-^ to 10-^ and 7 from 10^^ to 10. 

while in the second case it reaches outward and could be considered approximately constant 
within it. The Yukawa couphng was chosen so that the effect on Tc could be sizeable (that is, 
of order 0(10-^)). 

This enables us to build the exclusion plot of Figure 2, by imposing that AT^ < 4 x 10~^, 
the accepted bound derived from Solar luminosity constraints [13]. It is superimposed on the 
different bounds available [7]. 

3.2 VAMP models 

As previously discussed, we consider usual fermions with a variable mass term 5m given by 
the coupling to the quintessence-type scalar field; its mass is then given by the usual Higgs- 
mechanism related term m plus this new VAMP term, a small contribution: m — mniggs + ^fn- 
We adopt, however, a "worst case" scenario, in which the electron mass is not mainly due 
by the Higgs mechanism plus a minor VAMP sector contribution, but fully given by the said 
VAMP component alone. This implies that the cosmological expectation value should be weakly 
perturbed, so that the electron mass does not undergo large variations. 

It can be shown (see Appendix) that this variable term leads to a geodesic deviation equation 
of the form 
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(60 t.O"f^' ID'W 



Figure 2: Exclusion plot for the relative deviation from unperturbed central temperature Tg, 
for A ranging from 10~^ to 10~^, and 7 from 10~^ to 10 (tip at the top), superimposed on the 
available bounds [7] 



. !, . „ 2q; 

x^x" H x" 

a 



(26) 



Assuming isotropy, one obtains the acceleration 



m 



m _ 



<t>. 



a = ClNewton ~i fl'ftc^ X li^ l' = ^Newton ~i rQbc^ X Uj- -f . 

mm (p <p 



(27) 



where the prime and the dot denote derivatives with respect to r and t, respectively. Considering 
the Newtonian limit Qab — diag{l, — 1, — 1, — , 1) so that g^cX^x'^ = 1 — ~ 1, one finds a radial, 
anomalous acceleration plus a time-dependent drag force: 



— < 0, an 



(28) 



Notice that this radial Sun bound acceleration has the qualitative features for a possible anoma- 
lous acceleration measured by the Pioneer probes [9] . 

The time-dependent component should vary on cosmological timescales, and can thus can 
be absorbed in the usual Higgs mass term. Hence, one considers only the perturbation to the 
Lane-Emden equation given by the radial force, a a — 0'/^- 
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-GM{r) + qat' 



pdr 



(29) 



which translates into 



1 d 

J.2 fj^j. 



dP 
p dr 



-A-kGp - 



ct>{r) 



r 0(r) 



(30) 



Defining the dimensionless quantities 



U 



GM 



= 2.12 X 10"^ , 



(31) 



G-^ = (n + l)iVy("+^)wy("+^) , 



where 



one obtains the perturbed Lane-Emden equation: 



(32) 



,de 



-^"(0 - 



Gr, 



(33) 



u m 

The Klein-Gordon type equation for the scalar field, written in terms of the ^ variable is, 
inside the star, given by 



1 



no + ^0'(o 



-puo<i 



-(p+l) 



(0 + 



Apcg"(0 
1^ 



(34) 



where now the prime denotes derivation with respect to ^. The parameter p is the mean 
molecular weight and it is assumed that there is one electron per molecule, that is, that the 
star is composed by "Hydrogen" with a molecular weight p. 

Beyond the star, in the Klein-Gordon equation one has the coupling to the constant number 
density of fermions in the vacuum = = 3 m~^, 



1 



-puoc, 



(35) 



These equations constitute a set of coupled differential equations for 0(^) and 9(0 and, 
of course, the continuity of (f) across the surface of the star must be addressed. A complete 
derivation can be found in Ref. [5]. In here, a and pc depend on Wn, Nn and related quantities, 
which are evaluated after the solution 6{^) is known, together with M and R. Hence, one 
considers their unperturbed values for the Sun: a = R/^i^^ = 1.009 x 10® m and Pc = 1.622 x 
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10^ kg m~^. Moreover, current solar estimatives indicate that ~ 0.62 mp, being the 
Hydrogen atomic mass [13] . 

For simphcity, we deal only with the case of p = 1, as in Ref. [11]. Instead of the potential 
strength uq, we work with the potential energy density < 1 of the scalar field. Before 
presenting the obtained numerical solutions, we develop the expression for VLy = V{(f)c)/ Pcrit, 
with pcrit — 1-88 X 10~^^ h'^ g cm~^; in what follows we chose h = 0.71. Therefore 



Uq 



Uq 



Xn 



V 



+ Xn 



Xn 



V 



2^uoXnv , 



(36) 



which implies 



Uo 



AXn 



(37) 



V 



We now rescale the scalar field so to work with a dimensionless quantity $ = 0/0*, where 
0* is the cosmological vev obtained by assuming as reference values A = Qy — 1, 



0* Pcrit 

2nv 

Hence, the cosmological vev for general A and uq is given by 



(38) 



(f>c = 



^vp. 



crit 



2Xnv 



(39) 



so that $c = (pc/ (pc — I A. Thus, the Klein-Gordon equation has the following form: 



i) Inside the star 



^"(0 + ^$'(0 



Pcrit ^ 



Pcrit P' 



ii) In the vacuum 



S Pcrit'^ Pcrit 



The perturbed Lane-Emden equation assumes the form 



(40) 



(41) 



1 d 



-^"(0 - 



Cn 1 



2 $'2f£V 



(42) 



Since the perturbation on d{^) is to be shown to be small, one can take the unperturbed 
function 9(^) ~ 6*0(0 when solving Eqs. (40) and (41), and then introduce the obtained 
solution for the scalar field $(C) in Eq. (42). First, one assumes that the scalar field is given 
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by its cosmological vev perturbed by a small "astrophysical" , $a(0> contribution, $(^) 
Jly/A + $a(0- Hence, the Klein-Gordon equation becomes 



Pcrit A* 



Pcrit 



inside the star, and 



<(0 + -MO ^ 

S Pcrit 



Pcrit 



(43) 



(44) 



in the outer region. 

Substituting by the Sun values a = = 1.009 x 10^ m, Pc = 1.622 x 10^ kg m'^, 

11 0.62 rup, = 3 m~^, one gets 



KiO + faiO ^ 6.8A 



4.79 X lO^^e^iC) - A ( 1 - 



2A$.(0 ^ 



(45) 



inside the star, and 



<(0 + -^KiO ^ 6.8A 



A 1 



2A$,(0^ 



V 



(46) 



in the outer region. 

Numerical integration of these equations enables the computation of the central tempera- 
ture's relative deviation; as boundary conditions for $({) it is imposed that both the field and 
its derivative vanish beyond the Solar System (about 10^ AU). One can see by inspection that 
the solution $a(0 is practically the same, regardless of the value for Qy, as one always assumes 
$a(0 f^y/A. However, it is highly sensitive to A. 

Also, one must verify the validity of the condition ^ f^y/A for chosen Qy A values. 
For this, note that $o(0 evolves as A~^ and, as stated above, it is fairly independent of fly to 
a very good approximation. Hence, choosing a smaller value for fly amounts to reducing A, 
both by lowering the field $a(0 increasing its upper limit, Jly/A. 

By the same token, each value of fly corresponds a maximum allowed value for the coupling, 
Xmaxi^v)- One then uses these values to numerically obtain to first order solutions for 9{^) and 
0(^), for say fly — 0.1, 0.4 and 0.7, as presented on Figures 3 and 4. This enables one to extract 
the variation of the central temperature, Tc. The maximum allowed values for A (depending on 
the chosen fly) arc: for fly = 0.1, X^ax = 1-24 x 10"^^; for fly = 0.4, A < 2.45 x 10"^^ and for 
fly = 0.7, A < 3.3 X 10"^^ The limiting case fly = 1 yields A < 3.93 x lO"^"*. 

None of the presented curves exceed the maximum allowed variation for Tc of 0.4%: the 
maximum of STc — 2.82 x 10~^ occurs for fly — 0.7, A — 2.82 x 10~^^. Hence, the luminosity 
constraint is always respected and the bound one must respect is A < 10~^^. 
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e / e„ - 1 




Figure 4: Perturbed solutions for 9{^) ~ T(0, for the = 0.1, A = 1.06 x 10"^^ case; other 
solutions overlap. 
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3.3 The Pioneer anomaly 

Following the method outlined above, we look at the hydrostatic equilibrium equation for 
a constant perturbation: 



-GM{r) + qat' 



pdr 



(47) 



and therefore 



1 d 

J.2 ^j. 



r^dP 
p dr 



-A-kGp + 



2aA 



(48) 



The last term is a perturbation to the usual Lane-Emden equation, which is given by 



,de_ 



-r(0 + 



Qa 



27rGpca^ 

The factor in the perturbation term can be written as 



(49) 



1/2 (n+l)/2r, 



2TTGp^a = [{n + l)GK7rY'' p\ 
Following the same steps as before, and defining the dimensionless quantities 



(50) 



C-^ = ^(n + l)7riy„, (5 



^.^ = 3.65x10-3 a, , 

GM Qq 



(51) 



one obtains the perturbed Lane-Emden equation 



1 d 



,d9 
d^ 



-r(0 + /9<:^. 



(52) 



As previously, the boundary conditions for this modified Lane-Emden equation are unaf- 
fected by the perturbation: from the definition p = Pcd{C)i gets 6'(0) = 1; the hydrostatic 
equation Eq. (48), in the limit ^ ^ still imposes \d6/d^\^^Q = 0. In the present case one 
has only one model parameter, (3, which can be constrained by the same luminosity bounds as 
before. Solutions for this equation with (3 varying from IO'^^Qq to IO'^^Qq (the reported value 
is of magnitude op ~ 3 x lO^^^a©) enable one to compute the relative central temperature 
deviation as a function of /3, as presented in Figure 5. 

From Figure 5, one concludes that the relative deviation of the central temperature scales 
hnearly with a^, as 5Tc ~ aA/ao- Thus, the bound 5Tc < 4 x 10"^ is satisfied for values of this 
constant anomalous acceleration up to a Max ~ lO^^a©. The reported value is then well within 
the allowed region, and has a negligible impact on the astrophysics of the Sun. 
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Figure 5: Relative deviation from the unperturbed central temperature, for ga ranging from 
lO-^^a© to 10-^^00. 

4 Conclusions 

In this work we have studied solutions of the perturbed Lane-Emden equation for three diferent 
cases, related to relevant scalar field models. We obtain bounds on the parameter space of each 
model from Solar luminosity constraints. 

The exclusion plot obtained for a Yukawa perturbation produces no new exclusion region in 
the parameter space A-Xy- This results from the low accuracy to which the central temperature 
Tc is known, when compared to the sensibility of dedicated experiments [7] . Therefore, it is fair 
to expect that improvements in the knowledge of the Sun's central temperature could yield a 
new way of exploring the available range of parameters. 

For the VAMP case we have shown that the Yukawa coupling of the VAMP sector is con- 
strained to be A < 10~^^, and that the Solar luminosity constraint is always respected. The 
numerical analysis reveals that and A should satisfy the relation X/fly < 10~^^. 

It has also been shown that the scalar field acquires its "cosmological" value just outside 
the star, leading to no differential shifts of the particle masses in the vacuum; thus, there is 
no observable variation of fermionic masses, and hence no violation of the Weak Equivalence 
Principle. 

Finally, we found that an anomalous, constant acceleration such as the one reported on 
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the Pioneer 10/11 spacecrafts is allowed within the Sun for values up to IQ-^a©, thus clearly 
stating that the observed value op ~ lO^^^a© has negligible impact on the central temperature 
and other stellar parameters. 

5 Appendix 

In order to encompass models with a variable mass, we consider the generalized Lagrangian 
density 



(53) 



The function a is, in the homogeneous and time-independent case, identified with the square of 
the rest mass. We now deduce the Euler- Lagrange equation for the timelike geodesies. Notice 
that there is no right side terms because r is an affine parameter: 



= ^T7^ = arQnhX X + aQnh.rX X — [(X QacX j = 



= {(^fiQah + (^gah,c) x"'x^ - Sas'aci" " '^■ag ac,bX"' x^ - 2agacx"' 

= (—gab + gab,c) X'^x'' - —gacX" - ^g^cbX^X^ - 2gacX'' = 

\ a J a 

= {—gab + gab,c - '^gacb) X"-X^ - —gacX" - 2gacX°- = 

\ a J a 



gacx" + 



-( -u - \ - 

2 \gac,b ~r gbc,a gab,c) 2oi^^^ 



a 



x^x" + 2-gacX"' 
a 



T^g"^ {gac,b + gbc,a - gab,c) - TTd" gab 



x'^ + 



1 
2 

2a 



2a- 



x^'x' + 2-5^Vei^'' = 
a 



a 



.cd „ -a 



^ab n _ gad 



a 



x^'x' + — = 
a 



(54) 



In the isotropic, Newtonian case, one has a — Q;(r, r ~ i), and thus 



a! 



2a ^ 



a^aN + —gbcX x'-Ur v , 

2a a 



(55) 



where the prime denotes derivative with respect to the radial coordinate. 

Using gab = diag{l, — 1, — 1, — , 1) so that gbcX^x'^ = l—v^ ~ 1, one obtains a radial anomalous 
acceleration plus a time-dependent drag force: 



ar ^ 2a ^ 

OA = 77^ < 0, qd^ < 

2a a 



(56) 
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